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^ ■ ABSTRACT 

■ Context. Very old halo stars were previously found to show at least four different abundance ' anomalies' , which models of fast rotating 
QQ ' massive stars (spinstars) can successfully account for: rise of N/0 and C/0, low ^^C/^^C and a primary-like evolution of Be and B. 

! Here we show the impact of these same stars in the enrichment of Sr and Ba in the early Universe. 

Aims. We study if the s-process production of fast rotating massive stars can offer an explanation for the observed spread in [Sr/Ba] 

■ ratio in halo stars with metallicity [Fe/H]< -2.5. 

' Mettiods. By means of a chemical inhomogeneous model we compute the enrichment of Sr and Ba by massive stars in the Galactic 

Ol halo. Our model takes into account, for the first time, the contribution of spinstars. 

Results. Our model (combining an r-process contribution with a s-process from fast rotating massive stars) is able to reproduce for the 
^ ■ first time the observed scatter in the [Sr/Ba] ratio at [Fe/H]< -2.5. Toward higher metallicities, the stochasticity of the star formation 

O ^' fades away due to the increasing number of exploding and enriching stars, and as a consequence the predicted scatter decreases. 

I I Conclusions. Our scenario is again based on the existence of spinstars in the early Universe. Very old halo stars were previously 

O . found to show at least four other abundance "anomalies", which rotating models of massive stars can successfully account for Our 

1^ ■ results provide a 5th independent signature of the existence of fast rotating massive stars: an early enrichment of the Universe in 

C/^ ' s-process elements. 

■ 

■ Key words. Galaxy: evolution - Galaxy: halo - stars: abundances - stars: massive - stars: rotation - nuclear reactions, nucleosynthe- 
sis, abundances 

>■ 

. 1, Introduction strongly dependent on the properties of the first stars, such as 

^^f-^ ■ their mass spectrum and chemical composition. 

^ Soon after the Big Bang, the appearance of the first stellar ^-^^^^^^ formation of the first stars is technically very 

. generations drastically changed the course of the history of challenging, as one starts from large-scale cosmological simu- 

g ■ the Universe by ennchmg the pnmordial gas with elements j^^^^^^^ j^s ) j^at are then zoomed into the highest den- 

^ heavier than helium (referred to as metals) through both stel- ^< 10-4 ^^^^^ng a dynamical range of at least 

lar winds and supernova explosions. High-resolution hydrody- the resolution of these simulations increases, the con- 

^ namical simulations of the formatK)n of the first stars suggest ^^^^-^^^ ^^-^ ^^^^^ emerging) field change drastically. The 

^> these objects to have formed in dark matter mim-halos, and ^^^^ ^^^j^^^^ ^^^^ ^^^^^^^ j,^^ ^^^^ 

^ to have played a key role m the formatio n and properties of ^^^^ ^^^^ theoretical simulations predicted that pri- 

_ the first galaxies (see [ Karlsson et al. |m for a review). One ^^prdial sta rs wer e predominantly m assive feromm et al.lll999l 

C3 . of the key questions, drivmg enormous theoretical and obser- |Ahel et al.ll2()()2L iBromm & Larsonll2()()4l) . current results sug- 

- - ■ vational efforts, is whether or not primordial environment in ^^^-^ ^^^^ ^^^^ ^^^^^ ^^^^^^ ^^^^^ 

which the first stars were bom made them different from present- ^^^^^^^ ^^^^ 50 ^ ^ ^^^^ ^^^^ low -mass stars 

day stars. The answer to this question is crucial to understand ^^^^^ ^j^^ have formed in th e very early Univ erse (IStacv et al.l 

the impact of the first stars, both in terms of metal injection [2oTil[2oT2Llciark et al.ll20TTt[G^dfeLll20rni2^^ 

and energetic feedb ack, on the formation of the first galax- , ^ , . t • , i , . . 

ies (iBromm &Yoshidal 1201 II) . Major endeavors, such as the These state-of-the art simu ations have become sophisti- 

James Webb Space Telescope or the 40m European Extremely '^^t^'^' now able to resolve the accretion disk onto the 

Large Telescope, expected to detect the first galaxies, have been Proto-stdlar core. However, on one hand the complexities such 

pushing the frontiers of high-resolution hydro-dynamical sim- *e efi^ects of turbulence, magnetic fields and negative feed- 

ula tions with the aim to pr edict their luminosities and col- ^^^^ by the forming protostar HII region on the accretion disks, 

ors (iJoggerst&WhalenlM . These predictions, however, are *e other hand problems related to the enormous reso- 

lution needed to correctly follow the accr etion process, leads to 
a still controversial scenario (see Ferrara 20121 for a recent dis- 



emailto: cescutti@aip.de cussion). As a consequence, the current mass-spectrum of the 
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first stars is still not known. Very interestingly, the most recent 
simulations seem to converge to one result, namely: that there is 
sufficient angular momentum in the proto-cloud collaps e to yiel d 
rapid l y rota ting stars, near the break-up speeds (e.g. IStacy et al.l 
l2QTlll2Q12h . As we will see below, the completely independent 
approach of Galactic Archaeology, had already pointed to the 
same direction. 

It is now recognized that Galactic Archaeology off'ers a pow- 
erful alte rnative way of probing the nature of the first stellar gen- 
erations (iFreem an & Bland-Hawthomll2QQ2l; iBromm & Yoshidal 
[201 1: Tumlinson bOld) . thus providing invaluable constraints to 
the above-mentioned simulations. It is in the Milky Way and 
its satellites that the oldest and most metal-poor stars in the 
Universe are observable, born at times or equivalent redshifts 
still out of reach of even the deepest planned extragalactic sur- 
veys. 

Low-mass stars (with masses below 80% of a solar mass) 
have lifetimes comparable to the age of the Universe. We, there- 
fore, expect that in their atmospheres, the elemental abundances 
of the gas at the time of their birth are mostly preserved, off'ering 
a local benchmark to cosmology. These stars thus keep memory 
of the unique nucleosynthesis in the first stellar generations, pro- 
viding invaluable constraints on the masses of the first stars and 
therefore on their stellar yields - the very basic input for mod- 
eling subsequent stellar generations, leading to the abundance 
distribution prevailing at present time. 

Reliable age determination, especially for the oldest stars, 
has been one of the biggest challenges in astronomy, although 
a breakthrough is expected to come from astro- seismology 
(IChiappinill20 1 2lfl . In the meantime, the strategy adopted so far 
has been to look for the most metal-poor stars in the Milky Way 
and its satellites, presumably formed from a gas enriched only 
by a few first supernovae. 

The hunt for the most metal-poor stars in the Universe is not 
an easy task. This search is restricted to nearby (resolved) stel- 
lar populations: the MW halo and its satellites. The chemically 
most primitive stars are extremely rare. Despite their scarcity, 
currently, in the MW halo around 130 halo stars, with metallic- 
ities below [Fe/H] = -3, are known (Frebel 201 0, for a recent 
review). Only two of them h ave [Fe/H] < -5 ( Fr ebel et al.ll2005l 



l2008l:IChristlieb et al.ll2002 '): more recently, Caff'au et al.l (l201ll 



discovered a star with a [Fe/H] = -4.89 (so slightly above -5) 
but with the lowest global metallicity (Z< 7.4 x 10"^) being 
not enhanced in C and N diff'erently from the Frebel 's and the 
Christlieb's stars. 

Recently, another window into the very metal-poor uni- 
verse was provided by the discovery of the so-called ultra-faint 
dwarf spheroidal galaxies. These objects are dark matter dom- 
inated and have very low average metallicities (some show a 

mean [Fe/H] 2.6), thus off'ering an excellent environment 

in which to search for more fossil records of the early enrich- 
ment. The search for more of these Rosetta stones of the his- 
tory of the Universe constitutes the goal of several planned and 
ongoing surveys, for details the interested read er is referred t o 
SEGUE ( Yannv et al. 2009), APOGEE (Eisenst eliretaD|201lh . 
LAMOST (Zhao et al. 2006), SkyMapper ( Keller et al. 2007), 
4M0ST (de Jong et al. 2012), WEAVE (Balcells et al. 2010). 

The new information on cosmic abundances provided by the 
still modest samples of very metal-poor stars in the halo (still 
confined to mostly nearby and brighter objects) and the MW 



satellites have already proven to be a rich source of informa- 
tion on the nature of the fi r st stel lar generations. In their pio- 
neering work, Cay rel et aP (l2004h derived abundances of sev- 
eral elements for a fairly large sample of very-metal poor nor- 
mal starfl obtained with the VLT-first stars large-program. Their 
data of unprecedented quality revealed a striking homogeneity in 
the chemical properties of halo stars. 

At first sight, this uniformity was interpreted as a challenge 
to the view that the whole galactic halo formed from the suc- 
cessive swallowing of smaller stellar systems with independent 
evolutionary histories, as predicted from the cold dark matter 
theory. The situation turned out to be more complex, as the same 
stars w ere found to show la rge scatter in r- and s-p rocess ele- 
ments (iFrancois et al.l l2007l) and in the N/0 ratios (ISpite et al.l 
I2005h . 

The first stars are believed to have been radically diff'erent 
from present-day massive stars, as they were metal-free. The 
lack (or only traces) of metals leads to faster surface rotation 
velocities, as metal-poor stars are more compact than metal rich 
ones. Stars formed from a gas whose global metallicity is below 
-1/2000 that of the Sun (i.e. Z = 10"^) could attain rotational 
velocities of 500-800km/s (depending on the stellar mass, here- 
after called spinstars - see lMaeder & Meyneil2012l for a recent 
review). 

Two important eff'ects, of interest for the present work, arise 
in spinstars: a) rotation triggers mixing processes inside the star, 
which leads to the production of imp ortant quantities of primary 
^"^N, ^^C and ^^Ne ("Meynet et al." 2008, and references therein). 
The production of primary ^^Ne has an important impact on the 
s-process nucleosynthesis in spinstars compared to non-rotating 
stars, increasing by orders of ma gnitude t he s-process yields of 
heavy elements (e.g.lPignatari et al. 2008: IChiappini et 511201 ll: 
Frischkn echt et al ]'2012'); b) the strong mixing caused by ro- 
tation enriches the stellar surface, thus increasing the opacity 
of t he outer layer, and producing line driven winds (but see 
Muij res et al.ll2012l) . In addition, fast rotation can also lead to 
mechanical mass-loss. Both mechanisms are able to trigger non- 
negligible mass-loss already at very low metallicities (a result 
not achieved by standard models without rotation). 

We have shown that a generation of spinstars in the very 
early Universe is currently the most promising for the increas- 
ing ratio of N/0 ratio observed in very metal-poor normal stars 
in the Galactic halo (IChiappini et al.l l2QQ6k . Our models also 
naturally account for the observed incr ease in the C/0 ratio 
in halo stars towar ds low metallicities d Chiappini et al.l I2Q06I: 
Fabbian et al. 2009). Moreover, if spinstars were common phe- 
nomena in the early Universe, the metal-poor ISM of galaxies 
with a star formation history similar to the one inferred for our 
galactic halo, should have ^^C/^^C ratios between 30-300 below 
[Fe/H]=-3. Without fast rotators, the predicted ^^C/^^C ratios 
is 4,500 at [Fe/H]=-3 .5 increasing to 31,000 at [Fe/H] = -5.0 
(IChiappini et al .1120081) . Current observations of this isotopic ra- 
tio in very metal-poor giant stars in the galactic halo (ISpite et al.l 
l2006l) a gree better with chemical evolution models including fast 
rotators. Recently, the same results for C isotopes and for the 



^ In some rare cases it was possible to estimat e an age by measur - 
ing long-lived radioactive elements as Th and U dCavrel et al .1120011) . 
although with still non-negligible uncertainties. 



^ Around 20% of stars with [Fe/H] < -2.0 show extremely high 
C abundance, with [C/Fe] > 1.0. These stars are named C-Enhanced 
Metal-Poor Stars (CEMP). The remaining 80% stars are called normal 
stars (Lucatello et al. 2006). 

^ Indeed, asymptotic Giant Branch stars (AGBs), that produce ^"^N 
and ^^C efficiently, have not had time to contribute to the chemical en- 
richment of the early Universe due to their long lifetimes; Another vi- 
able explanation not fully explored yet is the contribution of SAGB. 
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ratio N/0 have been obtained in iKobayashi et alJ (1201 ll) using 
another chemical model. 

To test the latter predictions, challenging measurements of 
the ^^C/^-^C in more extremely metal-poor giants or tumoff 
stars are required (only feasible with 30m-class telescopes). 
Furthermore, we explain the apparently contradictory finding of 
a large scatter in C, N and the almost complete lack of scatter 
in [a/Fe] ratios of the same stars as being a consequence of the 
combination of a stochastic star formation rate in the early uni- 
verse (hence, incomplete sampling of the IMF and incomplete 
halo mixing) combined with a strong yield dependence on the 
stellar mass for these particular elements (Cescutti & Chiappini 
I2Q1 Q): spinstars provide a nucleo synthetic site in which the 
yields of C, N are strongly dependent on stellar mass. 

In summary, we claim the rise of N/0 and C/0 and the low 
^^C/^^C in very metal-poor halo stars to be three footprints of the 
existence of spinstars in the earliest phases of the Univer se. A 
forth signature was recently suggested by ( Prantzosll2012l) who 
argue that the observed primary-like evolution of Be and B can 
be explained if Galactic Cosmic Rays are accelerated from the 
wind material of rotating massive stars, rich in CNO, hit by the 
forward shock of the subsequent supernova explosions. These 
four footprints have been found in the very metal-poor Universe 
([Fe/H] < -3). 

In the present paper we suggest that two puzzling results in- 
volving n-capture elements such as Y, Sr and Ba might again 
point to the existence of spinstars in the early Universe, namely: 

- The large scatter in the abundance ratios of s-process el- 
ements from the first and second peaks, e.g. [Sr/Ba] (see 
Section O, which cannot be explained even by inhomoge- 
neous models that otherwise match the lack of scatter in a 
elements, and the large scatter of the same elements with re- 
spect to iron (e.g. Sr and Ba in Cescutti 2008); 

- The large enhancements of Ba, La, Y and Sr found in stars of 
NGC 6522, a bulge glo bular cluster shown to be the oldest 
one in the Milky Way (iBarbuv et al. 1 120091: IChiappini elaP 
I2OIII) . 

Indeed, IChiappini et al.l (1201 ll) have shown that: a) the scat- 
ter in the [Sr/Ba] and [Y/Ba] for very metal-poor halo stars is 
similar to that found in the NGC 6522 stars, and b) that the en- 
richment of Sr, Y, Ba and La in 5 out of 8 stars of NGC 6522 
could be due to either the s-process activation in early genera- 
tions of spinstars, or the s-process due to the contamination from 
a low-mass AGB stars in a binary system, with initial metallic- 
ity coincident with the metal content of the clustefl However, 
the remaining 3 stars showing the highest [Y/Ba] are not com- 
patible with s-process nucleosynthesis in AGB stars and can be 
readily explained from early enrichment from spinstars. It is 
now possible to investigate this hypothesis in a more quanti- 
tative way by using inhomogenous chemical evolution models 
computed with the most recent calculations carried out by our 
group (Frischknecht et al. 2012, Frischknecht et al. in prepara- 
tion). In the present work, we will focus on the observed scatter 
of [Sr/Fe], [Ba/Fe] and [Sr/Ba] in very metal-poor halo stars. We 
argue that fast rotating massive stars have also left a footprint in 
the early enrichment of heavy elements, such as Sr and Ba. In a 



In the first case the chemical enrichment is from stars polluting the 
primordial material before forming the cluster, whereas in the second 
case the chemical enrichment would happen via AGB -mass transfer at 
whatever point in the history of NGC 6522, thus polluting the stars we 
are now observing. 



forthcoming paper we will then specifically discuss models for 
the Galactic Bulge. 

The paper is organized as follows. In Section 2 we describe 
the current situation with respect to the predicted stellar yields 
of heavy elements, and justify the use of empirical yields for 
the r-process In Section 3 the observational data is described. 
Section 4 and 5 describe the chemical model and the adopted 
stellar yields, respectively. In Section 5 our results are presented 
and in Section 6 we summarize our conclusions. 

2. The production of heavy elements: an open 
debate 

Elements with Z >30 are labelled neutron capture elements: they 
are mainly formed through neutron captures, and not through 
fusion, because charged particle reactions on elements beyond 
iron (Z = 26) are endothermic 0. The neutron capture process 
is split in rapid process (r-process) or slow process (s-process) 
depending whether the timescale for neutron capture t„ is faster 
or slower compared to radioactive beta deca y r^, according to 
the initial definition by iBurbidge et al.l (Il957h . 

There are several sites of production for the s-process. 
Massive stars produce the so-called weak s-process (60<A<90) 
and intermediate and low mass stars produce the ma i n s-pr ocess 
(up to lead) - see recent review by iKappeler et al. l (l2Qnl) . For 
r-process, the final states of massive stars where extremely 
high neutron density are achievable, seem to provide a promis- 
ing site (although still very uncertain, see Thielemann et al.l 
,20111 A multiplici ty of hypoth esis do exist, such as n eutrino 
winds during SNII (IWooslev et al. 1994; Takahashi etal] |l994l: 
lThom pson2003: Arcones et al. 2007; Farouai et al. 2009), neu- 
tron stars mergers (Freiburghaus et al. 1999; Gorielv et al. 201 1; 
Korobkin et al. 2012), 0-Mg-Ne core exp losion (IWanaio et al. 
2011), asymmetric explosions (Cameron! 120031: lAmould et al. 
[2007), quark novae ( Jaikumar et al. 2007) and magnetorotation- 
ally driven supernovae ( Winteler et al. 2012*). 

The general complexity of these environments, coupled with 
the inaccessibility of the nuclear data of the isotopes involved, 
leads to still uncertain predictions for the stellar yields. Only re- 
cently a few SNII simulations have started to successfully ex- 
plode stellar cores and the results are still controversial. This 
is the case for th e simulations of Q-M g-Ne core collapse su- 
pernovae (SN) by iKitaura et al.l (|2006|) and the nucle osynthesis 
calculated on this basis bv'Wanaio e t al.1 (l201 lll2009l) . 

Even though the r-process is not clearly understood, most re- 
cent studies in this field point to massive stars (with the remark- 
able exception of the neutron stars mergers) as the most promis- 
ing production site. If this is the case, the chemical enrichment 
timescale of the r- and s- processes are distinct, namely: few 
tens of million years for the r-process nuclei (typical lifetime 
for a massive star), and more than 0.5 G yr for the bulk of pro- 
duction of s-process in AGB stars (K appeler et al.ll20T]1) . Some 
production of s-process is also expected to take place in massive 
stars, although the overall contribution is negligible, especially 
at the very metal-poor regime (see Raiteri et al. 1992). For this 
reason, it has been common to associate the n-capture elements 
in the extremely metal-poor (EM P) stars of the h alo to r-process 

production, as pointed out first by lTruran](ll98ll) . 

The Solar System isotopic abundances (see iGrevesse et al.l 

I2OIO ) can be split into r- and s-process contributions. 
Theoretically the s-process abundance pattern is known, and the 

^ we are omitting the p-process, responsible of less than 1% of heavy 
elements. 
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solar s-process pattern is obtained by scaling to the abundances 
of isotopes exclusively produced by s-process 0. Hence, the r- 
process contri bution to the solar system (Arlandini et al. 1999; 
ISimmerer et al.l |2QQ4|) . is determined by subtraction of the s- 
process contribution. In this way the so-called r-process con- 
tribution to the solar system is found, although not free from 
theoretical uncertainties in the s-process contribution, chemical 
evolution, nor other processes that could be hidden. 

Although the relative amount of neutron capture elements is 
small compared to those of the lighter elements (roughly 10"^), 
their absorption lines in the stellar spectra can be observed and 
their abundances c alculated even for th e most met al-poor stars. 
Since the w ork by ISneden et al.l (Il996l) (and also ISneden et al.l 
l2QQ3ll20Q8l) a remarkable agreement between the abundance pat- 
tern from B a up to Th measured in a very r-process rich EMP 
star (CS 22892-052) and the r-process contribution to the solar 
system has been found. This appeared to be a confirmation for 
the theory of the s-process from one side and of the idea that r- 
process was the only producer of (heavy) neutron capture in the 
early universe. 

On the other hand, the same invest igators ISneden et al.l 
(Il996h (and also ISneden et al.ll2003l l2QQ8l) have shown that the 
pattern for light neutron capture elements (from Sr to Ag) not 
always matched the so-called r-process solar pattern. Moreover, 
the ratio of light to heavy neutron capture elements in other EMP 
stars were not constant, as expected for two species coming from 
the same process. The light elements are more enriched in some 
EMP compared to the expectations of a pure r-process pattern, 
and show an increasing scatter towards low metallici ties, up to 
[Fe/H] ~ -4 (iHonda et al.ll2004al:lFrancois et al.ll2007l) . 

A new process was invoked to provide the additional contri- 
bution for light neutron capture elements and, at the same time, 
possibly produce the spread in the ratio of light to heavy neu- 
tron capture elements. Travaglio et al. (2004) called this process 
lighter elements primary process (LEPP), without specifying a 
possible r- or s- origin (although typically s-process are not a 
primary processfl Later jMontes et al.l (l2007l ), analyzing the pat- 
tern of this possible missing process, obtained two viable solu- 
tions, one with typical neutron fluxes of r-process, and the other 
closer to s-process values, both of them able to provide roughly 
the correct pattern. 

For the first time in the work by iPignatari et al.l (l2008l) . it 
was shown, assuming a primary amount of ^^Ne guided by rotat- 
ing models, that rotating massive stars could produce significant 
amount of s process at low Z. More recently. iFrischknecht et al.l 
(12012) calculated full stellar evolution models using an extended 
nuclear reaction network and showed that fast rotation can boost 
the production of s-process at low metallicities. 

Here our goal is to test the impact of these yields in the 
early chemical enrichment of the Galaxy. Clearly, the s-process 
in fast rotating massive stars will not be the only contributor to 
the Ba and Sr in the early Universe, as most of it shall come from 
the still uncertain r-process. This means that in our models, we 
need to also include a production via r-process. Unfortunately, 
given the uncertainties discussed before, only few groups do pro- 
vide r-process stellar yields, namely: for neutrino wind models 
(Arcones & Mon tes 20 U ) and for 0-Ne-Mg core collapse mod- 
els ( Wanajo et al.ll201 ll) 7However, as recently shown by Hansen 
et al. (2013) for the case of Sr, these predictions can still vary by 



^ Isotopes shielded against r-process by the corresponding stable iso- 
bar of Z-2 

^ Note that IWasserburg et al.l (Il996l) had already invoked two r- 
processes to explain radioactive isotopes anomalies in solar meteorites. 



more than 3 dex, and although the models quoted above do pro- 
duce Sr, they do not produce Ba. Other processes, more promis- 
ing for Ba production, such as neutron star merging or magneto- 
rotationally driven SN have not yet led to detailed stellar yields 
computations. 

In the present work we follow an alternative way to deal with 
the problem of the absence of robust theoretical yield predic- 
ti ons for the r-pro cess, already adopted by other au thors such 
aslTravaglio et al.l (Il999ll2004 . llshimaru & Wanaiol (Il999h and 
ICescutti et al.l (20oK namely: we adopt empirical yields, which 
are extracted from the available observational constraints (see 
Section 4.2.1). We notice that it turns out that the empirical 
yields we estimate are within the uncertainties of the present r- 
process models. This will enable us to quantitatively estimate the 
impact of spinstars on the early enrichment of heavy elements, 
as well as guiding theoretical eff'orts on r-process nucleosynthe- 
sis. 



3. Observational data 

3.1. Metallicity distribution function 

The chemical evolution model assumptions on star formation 
history and winds are such that the overall shape of the metallic- 
ity distribution (MDF) of the halo stars is reproduced. The MDF 
is strongly aff'ected by the selection criteria of each sample, and 
most studies have t ried to take this in to account. A comparison 
with the results by lYong et all (|2012|) . who focused on the ex- 
tremely metal-poor tail of the M DF shows a good agreemen t 
with the theoretical predictions of ICescutti & Chiappinil (l2010l) . 
which is essentia lly the same model adopted here (see Fig. 3 of 
lYong et al.ll20T2b . 

Two other recent analysis making use of data from the 
Hamburg/ESO objective-prism survey (HES) (Christlieb e t al.l 
2008) have been recently carried out wit h the goal to d etermine 
the whole halo MDF, namely: the Schorck et al.l (12009.) MDF 
based on metal-poor giants and the MDF from iLi et al.l (1201 Oh 
based on main sequence turn off' stars (MSTO) stars. Both anal- 
ysis use statically well understood selection criteria to obtain a 
large number of metal-poor giants (1638) and metal-poor MSTO 
stars (617) respectively. 

Here we compare our model predictions with the MDF 
traced b y MSTO stars rather then giants because, as shown by 
|Li et al.l (2010), the dwarfs are less aff'ected by survey-volume 
corrections. We are aware that this sample lacks the extremely 
metal-poor tail of stars with [Fe/H] < -3.5, but larger statisti- 
cally complete samples are required for this purpose; fortunately, 
such samples will be obtained from much larger and deeper 
surveys in the near future, such as from SEGUE-2, LAMOST, 
SkyMapper and 4M0ST,. 

Concerning the interesting results by'Carol lo et al.l(l2007l) re- 
garding the possible duality of the halo, in the present work we 
concentrate only on the chemical evolution of the inner halo to 
which our data belong. 



3.2. Chemical abundances 

We adopt observational abundance ratios from literature; the 
data for the neutron capture elements and for the (^-elements are 
those collected by Frebel ( 2010)0, labelled as halo starfl. 



^ http://cdsarc.u-strasbg.fr/cgi-bin/qcat7J/AN/33 1/474 
^ The list of the a ut hors we u s e from the collecti o n are : 
iMcWilliam et"aD ([T995h : iMcWilliamI ([19981) : IWestin et al.l ([200Qh : 



4 



Cescutti et aL: The 5th signature of spinstars in the Galactic halo 



Among the halo stars collected, we differentiate the normal 
stars from the carbon enriched metal- poor (C EMP) stars. We 
follow the definition given by Masse ron et alJ (^2010) where a 
CEMP star is defined as having [C/Fe]>Q.9 . The latter can be 
sub-classified as (see iMasseron et al]|2QlQl) , based on Ba and 
Eu, namely: 

- CEMP-s: [Ba/Fe]>l and [Ba/Eu] >0 

- CEMP-rs [Eu/Fe]>l and [Ba/Eu] >0 

- CEMP-no [Ba/Fe]<l and no Eu value 

- CEMP-r [Ba/Fe]<l and [Ba/Eu] <0 

- CEMP-low s [Ba/Fe]<l and [Ba/Eu] >0 

We note that in our sample we have just one CEMP with 
[Ba/Fe]<l with Eu measurement, which is classified as CEMP-r 
and so the CEMP-low s category is missing. For our purposes, 
the above categories are not strictly necessary, but we prefer to 
refer to literature than create "ad hoc" definitions. The impor- 
tant distinction in the present work is between CEMP-(r)s and 
CEMP-no, so whether a strong signature of s-process is present 
or not. 

Indeed CEMP-s most likely originat e from binary mass 
transfer froni a p revious AGB companion (iBisterzo et al 11201 21: 
iLugaroetaT 2012) and for this reason CEMP-s do not reflect the 
chemical evolution and hence they are not reflecting the chemi- 
cal evolution of the ISM. 



4. The chemical evolution models 

Homogeneous chemical evolution models alone cannot take ad- 
vantage of the precious information recorded in the scatter of 
key element abundance ratio of old stars, in difl'erent environ- 
ments. The need for inhomogenous chemical evolution models 
(i.e. models that relax the instantaneous mixing approximation), 
in turn, requires a cosmological framework, in which it is pos- 
sible to predict the properties (volumes) of the inhomogeneities 
arising from both merger building blocks and unmixed ISM (e.g. 
feedback). However, in this case other uncertainties related to the 
field of cosmological simulations will play a crucial role in the 
results. Although we have started to work in this direction, we 
still think that it is important to develop more simple inhomoge- 
nous models, of the kind of the one presented here. 

Indeed, these models will serve as a test-bench to study the 
diff'erent nucleo synthetic prescriptions proposed by the diff'erent 
scenarios. The goal is to identify key abundance ratios, as well 
as the most promising stellar yields, to be implemented in cos- 
mological simulations. This step is crucial to avoid spending un- 
necessary CPU time with models that do not meet the minimum 
requirement set by a first comparison with observations. 

The neutron capture elements spread in EMP stars was 
already the su bject of other investigations. We recall here 
the w o rks of iTsui imoto et al. ( 1999), Ishimaru & Wanajol 
il999'), "TravagHo et al. (2001), Argastet al., (i2002l 2004) and 
iKarlsson & Gustafsson ( 2005). Most of these models repro- 
duced the spread in n-capture elements observed in EMP stars, 
but left unexplained other problems: a much smaller scatter 
observed for the other elements (for in stance a-elements, see 
ICavrel et al.l 120041: iBonifacio et"al]l2012D . the reduction of the 

'Aoki et al. (2002); Cowan et al. (2002); Ivans et al. (2003); Honda et alj 
(2004b); Aokietal. (2005); Barklemetal. (2005); Aoki et al. (2006); 
Ivans etal. (2006); Masseron et al. (2006); Preston et al. (2006); 
Aoki et al. (2007); Francois e t al. (2007):,Lai et al. (20()7): Cohen et al 
(2008); Lai etal. (2008); IRoederer et all (l2008h : IBonifacio et all 
(>2009i) : .Havek et al.. aOO% 



spread toward higher metallicity for the neutron capture ele- 
ments and the large scatter observed in N/0 and C/0. 

This has pointed to a diff'erent solution in which the inhomo- 
geneities arise from a stochastic formation of massive stars. The 
spread can be generated by the enrichment of diff'erent species 
if they are produced by diff'erent ranges of masses, as show n for 
heavy neutron capture elements as a function of iron in .Cescuttil 
( 2008); at the same time the spread is reduced after few genera- 
tions of stars and it is small for elements which share the same 
range of production (as a-elements and ir on). This approach has 
also been used for the CNO elements (see ICescutti & Chiappinil 
2010) to investigate the impact of the inhomogenous modeling 
on the observed spread (in particular of the ratio N/0) in EMP 
stars. 

Here we consider the same chemical evolution model 
adopted in Cescutti & Chiappini (2010), based on the inhomoge- 
nous model developed bv ICesciittil (I2008h and on the homoge- 
neous model of Chiappini et al.l (l2008l) . 

We consider that the halo consists of many independent re- 
gions, each with the same typical volume, and each region does 
not interact with the others. We decided to have a typical vol- 
ume with a radius of roughly 90 pc and the number of assumed 
volumes is 100 to ensure good statistical results. The dimension 
of the volume is large enough to allow us to neglect the interac- 
tions among diff'erent volumes, at least as a first approximation. 
In fact, for typical ISM densities, a supernova remnant becomes 
indistinguishable from t he ISM - i.e., merge s with the ISM - 
before reaching ~ 50pc dThomton et al.|[l998l) . less than the size 
of our region. We do not use larger volumes because we would 
lose the stochasticity we are looking for; in fact, larger volumes 
produce more homogeneous results. 

In each region, we assume the same law for the infall of the 
gas with primordial composition, following the homogeneous 
model by Chiappini et al.. (.20081 so we have the following gaus- 
sian law: 

dGaSinit) ^ ^-it-t^fial 
dt 

where to is set to lOOMyr and (Tq is 50Myr. Similarly, the star 
formation rate (SFR) is defined: 

SFR{t)^{pgas{t)t' (2) 

where PgasO) is the density of the gas mass inside the considered 
volume. Moreover, the model takes into account an outflow from 
the system: 

oc S FR{t) (3) 

dt 

So, knowing the mass that is transformed into stars in a timestep 
(hereafter, M"fJ^), we assign the mass to one star with a random 
function, weighted according to the initial mass function (IMF) 
of ScaLa (1986) in the range between 0.1 and lOOM©. Then we 
extract the mass of another star and we repeat this cycle until 
the total mass of newly-formed stars exceeds M^fJJ!^. In this way, 
in each region, at each timestep, the MJfJJ!^ is the same but the 
total number and mass distribution of the stars are difl'erent, and 
we know the mass of each star contained in each region, when 
it is born and when it wil l die, assuming the stellar lifetimes of 
"Mae der & Me vnet* ('1989'). At the end of its lifetime each star 
enriches the ISM taking into account the enrichment due to that 
star and due to the already present mass of each element locked 
in that star when it is bom. 
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Fig. 2. [Ca/Fe] and [Si/Fe] vs [Fe/H], from left to right. The density plot is the distribution of simulated long living stars for our model; 
the density is in a logarithmic scale, normalized to the peak of the distribution and the bar over the plot describes the color scale assumed. 
Superimposed to the density plot, we show the abundances ratios of halo stars (data from .Frebel^O 1 0.) . 




Fig. 1. Comparison between the model MDF in blue line and the 
observed halo MDF by Li et al. (2010): main-sequence turnoff 
stars in the HESS (Hamburg ESO) for which [Fe/H]< -1. 



Yields for a e lements and Fe are the same as in 
iFrancois et al.l (|2QQ4|) . The model takes into account the produc- 
tion by s-process from 1.5 to 3Mo stars and SNIa enrichment. 



as inll Cescutti et al.l (l2Q06l) . The model does not include the most 
recent and accurate yields for AGB, and the very low metallic - 
ity yields for AGB are extrapolated. However, due to the longer 
lifetime compared to the halo star formation history, this produc- 
tions is hardly seen in the results of the model for [Fe/H] < - 1 
as well as in the abundance pattern of normal EMP. Super- AGB 
stars can be also involved in the production of neutron capture 
elements but at the present time yields for those elements from 
SAGB are not available; however, the production of s-proces s el- 
ements in their convective pulses is expected to be small ( SiessI 
[2010) . Additional studies are required to clarify this aspect. 



4.1. First test: reproducing the MDF and the [a/Fe] low 
scatter 



This m odel is able to w ell reproduces the MDF measured for the 
halo bv lLi etal.l(|2QTQl) based on HES data of main sequence turn 
off stars (see Fig.[T]). This comparison shows that the timescale 
of enrichment of the model looks like that of the halo stars in the 
solar vicinity. 



In Fig. [2l we show the predictions of our model for the a 
elements Ca and Si. Our model predicts a small spread for these 
abundance ratios, compatible with the observational data. In ad- 
dition, our model predictions show a slightly increase in the scat- 
ter in the very-metal-poor regime. Interestingly, inhomogenous 
models, as the one shown here, do predict the existence of a few 
outliers with low values of [a/Fe] . 



6 



Cescutti et aL: The 5th signature of spinstars in the Galactic halo 




[Fe/H] 

Fig. 3. [Ba/Fe] vs [Fe/H] abundances ratios of the halo stars (data from ^ Frebelll20 1 Ol) : black open circle are normal stars, blue open 
circle for CEMP-s stars (open triangles without Eu value) magenta open circle for CEMP-rs, cyan filled circle for CEMP-no and 
red X marker for the CEMP-r star. The error bars represent the mean and the standard deviations for the normal stars abundances 
calculated over diff'erent bins in [Fe/H] . The bins are calculated in such a way that each bin contains the same number of data. The 
results of the homogenous model with the assumed empirical yields for Ba is shown by the solid line; the dashed line shows the 
results when only the "standard" r-process site is considered. 




Fig. 4. Yields used, normalized to the solar (Grevesse et al. 2010), for the ratios of [Ba/Fe] and [Sr/Fe], as a function of the stellar 
mass and the metallicities 



7 



Cescutti et aL: The 5th signature of spinstars in the Galactic halo 



4.2. Stellar yields for heavy elements 

4.2.1 . Empirical yields for the r-process 

As mentioned in Section [2l the site of production of r-process 
elements is still a matter of debate. In particular, most of the 
neutron reactions involving r-process elements are out of reach 
of nuclear physics experiment. In such a situation, we must use 
observational data to guide the theoretical models for this pro- 
cess. 

Given the ab ove we adopt here the following approach (see 
also ICescutti et al. 2006), namely: we compute a homogeneous 
chemical evolution model where the yields of Ba are chosen 
such as to reproduce the mean trend of [Ba/Fe] versus [Fe/H] 
(see Fig. [3]). The resulting yields are what we call here empirical 
yields. 

The empirical yields are obtained as the simplest array able 
to reproduce the observed trend of increasing [Ba/Fe] with 
metallicity, but also present one important property, namely, the 
need for two different sites of production (see Fig. I?]): 

- a strong production in a narrow mass range S-IOM© that we 
call standard r-process site; 

- an extended contribution coming from stars with larger 
masses (from 10 up to SOM©), whose contribution is lower 
by ~ 2 dex than that of the standard r-process, and which we 
name here as extended r-process site. 

The standard r-process assumptions fit well the general con- 
cept of r-process production, so an intense production in relative 
rare events. On the other hand, the extended r-process seemed 
a necessary step to explain the presence of stars with [Ba/Fe] 

0.7 at metallicities lower than [Fe/H] 3.5 (Francois et al. 

2007, see Fig.©. This effect can be better appreciated from the 
inhomogenous results and it is a viable way to explain the ab- 
sence of stars lacking Sr and Ba (lRoedererll2013i) as in the case 
without the extended r-process). 

We here assume the r-process yields of Sr scale with those of 
Ba. The scaling factor adopted (of 3.16) has been taken from the 
solar system r-p r ocess contribution (see Sect|2) as determined by 
ISimmerer et al.l (12004 ^. Another meaningful assumption could 
have been the choice of the mean observed ratio between Sr and 
Ba in r-process rich stars (see Cowan et al. 2006). In this case the 
Sr/Ba ratio would be around 0.2 dex higher. On the other hand, 
this same ratio measured in the r-process rich star CS22892-052 
is close to the value that we have adopted. We stress that we have 
here selected a simple approach to the r-process production, al- 
though many authors have recently started to s uspect that this 
assumption may be not completely valid (Roederer et al.ll20TQt 
iBoyd et aT]|2012h : on the other hand, our simple approach is rea- 
sonable, since the emphasis of this paper is not on the r-process 
issue but on the contribution from spinstars to the production of 
strontium and barium. 

Interestingly the stellar mass range for standard r-process is 
close to the one predicted by theoretical models of ONeMg core 
supernovae (with initial masses toward the lower end of th e mas- 
sive s tars). Notice, however that the latest models (Wanaio 'et al.l 
1201 IL [20091) do not succeed in producing heavy neutron capture 
elements such as Ba. It is also worth mentioning that the stellar 
mass range 8< M© < 10 is in the transition between the super 
AGB stars and the electron capture supernova. For this reason 
we stress that the f inal fate for these stars is still rather uncertain 
(ISiessll2007Ll2QTQl) . 

(see comparison of model results with and without the extended r- 
process in Fig.|7]in the appendix) 



As we will see, the role of the extended r-process contribu- 
tion, that we need to take into account to explain empirically the 
observed [Ba/Fe] and [Sr/Fe] in EMP halo stars, can be replaced 
by the contribution of spinstars. 

4.2.2. The contribution of spinstars 

To illustrate the impact of rotating stellar models to the chemical 
enrichment of Sr and Ba in the earliest phases of the Universe, 
we now focus on three sets of inhomogeneous chemical evolu- 
tion models, computed with the following sets of stellar yields 
(see Fig. [4] and Table 1): 

- r-model: assume only empirical yields for the r-process 
(standard -\- extended); 

- as-model: assume only the standard r-process yields (no 
extended r-process contribution), plus the s-process yields 
coming from rotating stellar models (see below) 

- fs-model: similar to the model above, but with s-process 
yields coming from s-process in fast rotating stellar (spin- 
stars) models. 

The nucleosynthesis adopted in the as-model for the s- 
process come from unpublished results by Frischknecht (201 1, 
PhD thesis). In this set of yields, the s-process for massive stars 
is computed for Vm//V cnY=Q-4 and for a standard ch oice for the re- 
action ^^0(a, y) from lCaughlan & Fowlej (Il988h : they are com- 
posed of a grid of 4 stellar masses (15, 20, 25 and 4OM0) and 3 
metallicities (solar metallicity, 10"^, 10"^); in Fig.|4]we show the 
yields for the 2 lowest metallicity cases. We do not extrapolate 
the production toward stars more massive than 4OM0 (although 
it is realistic to have a production also in this range), but we ex- 
trapolated the Z = 10"^ grid down to Z=0. In addition to the 
s-process in massive stars, we take into account our empirical r- 
process enrichment but coming only from the standard r-process 
site. In this way we are decoupling the sites of production for the 
two processes. Interestingly, this figure shows that the s-process 
of spinstars plays here the role of the "extended r-process" which 
is not used in this model. 

It is worth pointing out that the yields we have used for 
the as- model are very con s ervati ve, among the models com- 
puted by lFrischknecht et aP (1201 2h for a 25Mo. In the fs-model, 
we have adopted a more robust production of s-process, which 
is achievable if we consider the s-process yields of massive 
stars computed for v/„//VcnY=0.5 (fast rotators) and for a reac- 
tion rate ^^0(a, y)^^Ne one tenth of the standard choice. Indeed, 
the adopted v/„//VcnY=0.4 in the as-model is lower than the 
value we have used in previous studies of spinstars (0.6-0.8) 
(Chiappini et al. 2006, 2008). I n ad dition, recent star formation 
simulations ( Stacy et al. 20 lH l2012i) support high rotation ve- 
locities compatible with larger Vini/Vcrit ratios. Further more, the 
current uncertainties in the ^^0(a,y) rate are still large. The 
adopted value of one tenth of the standard one is well within 
both the theoretical and recent experiments uncertainties. 

We do not have a fully computed grid for these parame- 
ters, but we have scaled the previous yields guided by the re- 
sults obtained w ith these parameters for a 25 M© of Z = 10"^ by 
iFrischknecht et al.l (120121) (cfr. in their paper Table 2), and ap- 
plying the same factor to all the masses at Z = 10"^, as shown in 
Fig. m Again, for this fs-model, in addition to the s-process in 
spinstars, we take into account the contribution by our empirical 
standard r-procQ^ss. 

Finally we note that IFrischknecht et al.l (1201 2l) provide only 
the pre- supernovae yields. However, the supernova shock only 
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affects the bottom of the carbon shell, and not the whole s- 
process rich regioiC3- As a consequence the bulk of the s-process 
rich material remains almost unaffected by the SN explosion (as 
is the also the case of light elements such as C, N and - see 
IWooslev et aDl2QQ2h . A more recent study bv iTur et alJ (l2QQ9h 
shows that certain isotopes can be affected by both explosive nu- 
cleosynthesis and also shell mergers possibly taking place dur- 
ing the early collapse. These changes will, however, not affect 
elemental ratios like [Sr/Ba] and the conclusions of our paper. 

5. Results 

Our main results are summarized In Figj5] We start by analyzing 
the results of the r-model (upper panel), which assumes only 
the contribution from massive stars via our empirical r-process 
yields. The spread obtained by this model matches the dispersion 
of [Ba/Fe] (upper, left panel), confirming that the hypothesis of 
a contribution by two distinct sites seems to be a good one to 
explain the data. In the particular case of this model, the two- 
sites of production are illustrated by a standard r-process which 
takes place in the lower mass range of the massive stars, and an 
extended r-process taking place in more massive stars. The first 
one is assumed to be much more efficient (larger quantities of 
ejected material) than the second. 

We consider also the match with the observed [Sr/Fe] ob- 
tained by the same model a good match, especially at very low 
metallicities: more than 80% of the stars below [Fe/H]< -3 are 
within the model results. As the Sr yields adopted here are ob- 
tained just by usi ng the Ba/Sr ratio ni atching the observed solar 
system r-process (ISneden et al.ll2QQ8l) . the results for this ratio 
are simply constant with metallicity (as shown in the upper, right 
panel). This suggests that some other physical process, taking 
place in the same mass range of what we have called extended 
r-process might be contributing to produce part of the Sr and Ba 
in the very early Universe (as the scatter is larger between metal- 
licities -2.5 and -3.5). In other words, the r-model is useful to 
highlight the issue we are trying to solve in the present work. 

Let us now turn to the results obtained with our as-model 
(see Table 1 and Fig. [3 second row), and see the impact of the s- 
process production of Ba and Sr by rotating massive stars. In this 
case we added a rather conservative s-process production from 
rotating stars, and turned off the contribution of what we have 
called extended r-process. The results for [Ba/Fe] and [Sr/Fe] 
are not completely satisfactory as the model cannot reproduce 
the low [Ba/Fe] ratios observed in extremely metal-poor stars 
with [Fe/H]< -3 (left and middle panels of the second row). 
This happens because this rather conservative model predicts not 
enough amounts of Sr, and particularly, of Ba (this can be seen in 
Fig. m by comparing the stellar yields of the extended r-process 
- dotted line - with the yields of the as-model at Z = 10"^ - 
dashed line). 

Despite the shortcomings described above, the interesting re- 
sult of this model resides in the [Sr/Ba] plot (second row, right 
panel), where it is clear that this new process produces an over- 
abundance and create an spread in [Sr/Ba] at [Fe/H] 3. From 

this one can conclude that the s-process produced in fast rotating 
massive stars seems to act in the correct direction. Nevertheless, 



The s-process takes place at the end of core He-burning and the start 
of shell carbon burning. At low Z, as explained in Frischknecht et al 
2012, the contribution from carbon shell burning is very small because 
of the strong neutron poisons during this phase. The composition of the 
s process rich layers are thus set by the end of He-burning and very little 
happens during the advanced stages at very low Z. 



these results are affected by the ratio of the yields of Sr and 
Ba predicted by this specific stellar evolution model (similar 
to model Bl of Frischknecht et al. 2012). These prescriptions 
produce [Sr/Ba] >2, whereas the EMP stars show at maximum 
[Sr/Ba]~1.5. 

Finally for the fs-model (Fig. [3 bottom panels), where 
we adopted less conservative stellar models predicting even 
larger s-process enrichment (essentially due to a larger Vini/Vcrit 
and a lower ^^0(a,y) reaction rate simil ar to the model B4 
of Frischknecht et al. 2012 - see Section l42l) , the theoretical 
predictions and observations show an striking agreement. This 
model not only reproduces well the [Ba/Fe] and [Sr/Fe] scatter 
(left and middle bottom panels), but can also account for the ob- 
served [Sr/Ba] spread (right bottom panel), at the right metallic- 
ity interval. This shows that with a less conservative production 
of s-process in fast rotating massive stars (as is the case in the 
fs-model), this process could play the same role as our extended 
r-process in the r-model. In addition, it shows that a standard 
r-process (taking place in the lower mass range of the massive 
stars) enters at play with a weight that increases as the metallic- 
ity increases (still in the very metal-poor range). 

In the scenario of the fs-model, the stars with high [Sr/Ba] 
abundances should show the pollution by the s-process in spin- 
stars. Hence, these stars should present a high value of the 
even/odd isotopes of Ba (as it is typically for an s-process). 
Interestingly, the pione ering measurements of Ba isotopes in six 
EMP stars by Gallaghe retaD (l2010i l20T2h have shown exactly 
the feature we predict with the above models, namely, a high 
level of even isotopes, compatible with an s-process production. 
It is important to note that four of these stars, with Sr abundance 
available in literature, present high [Sr/Ba] ratios, excluding an 
AGB pollution; for the remaning two stars a similar conclusion 
is supported by the Y abundance. We are aware that even with 
the highest resolution and signal to noise available, the measure- 
ment of the ratio of odd to even isotopes of Ba could still suffer 
from large uncertainties (such as broadening of the studied lines 
due to the differential hyperfine splitting of its isotopes). More 
measurements would be needed to confirm our suggestions. We 
stress, however, that at present, only our scenario can explain 
such observations. 

In Fig.O w e plot the [Sr/Ba] vs [ Ba/F e] similarly to the plot s 
in the works bv iMontes elaP (l2007l) and lFrancois et al.l (l2007l) . 
Our model can reproduce the peculiar behavior seen in EMP 
stars, namely, a high ratio of [Sr/Ba] at low [Ba/Fe], and the 
correct amount of scatter. This behavior in the model is deter- 
mined by the fact that by including the s-process production in 
spinstars we are able to decouple the sites responsible for the 
r-process process enrichment (S-IOM©), from the producers of 
s-process in the early Universe (i.e. the most massive spinstars, 
which are responsible for the largest contributions). In volumes 
where the most massive stars contribute more, the production of 
s-process will be important and the [Sr/Ba] ratios will be high. 
However, in this case the net enrichment in Ba and Sr will be 
low (due to the low yields produced by this process), which will 
then lead to low [Ba/Fe] and [Sr/Fe] ratios. On the other hand, 
in volumes where the contribution of the less massive stars is 
larger, the standard r-process will tend to produce larger levels 
of [Ba/Fe] and [Sr/Fe] enrichment, whereas the [Ba/Sr] will tend 
to be lower. The different mixture of these two limits, is dictated 
by the stochasticity of the star formation rate of each indepen- 
dent volume, which then lead to the observed spread. 
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Table 1. Nucleosynthesis prescriptions for the three studied cases. 



Model name 


panels in Fig.|5j 


s-process 


r-process 


r- 


Upper 


No s-process from massive stars 


standard + extended r-process site (8 - 30 M©) 


as- 


middle 


average rotators {Vinilvcmc - 0-4) 


standard r-process site (8-10 M©) 


fs- 


lower 


fast rotators {VinilVcmc - 0-5) 
and 1/10 for 0{a,y) reaction rate 


standard r-process site (8-10 M©) 




[Fe/H] [Fe/H] [Fe/H] 



Fig. 5. On the left [Ba/Fe], on the center [Sr/Fe], on the right [Sr/Ba], vs [Fe/H]. Upper, center and lower panels for "r-","as-" and 
"fs-"model respectively. The density plot is the distribution of simulated long living stars for our models; the density is in a log scale, 
normalized to the peak of t he distributio n (see bar over the Fig. [21 for the color scale). Superimposed, we show the abundances ratios 
of the halo stars (data from lFrebelll20 1 Ol) . The symbols are the same as in Fig. [3] 
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GO 




Fig. 6. [Sr/Ba] vs [Ba/Fe], the density plot is the distribution of 
simulated long living stars for fs-model; the density is in a log 
scale, normalized to the peak of the distribution (see bar over the 
Fig. [2] for the color scale), superimposed to the density plot, we 
show the abundances ratios o nly of normal halo stars (without 
carbon enrichment), data from iFrebell (|20 1 Ol) . 



6. Discussions and conclusions 

We have developed an inhomogeneous chemical evolution 
model for the halo with the aim of explaining the observed scat- 
ter (or lack of) in the abundance ratios of key chemical ele- 
ments in very metal-poor stars. The models presented here serve 
as a test-bench to study the different nucleo synthetic prescrip- 
tions proposed by different stellar evolution groups. The goal is 
to identify the most promising stellar yields for key abundance 
ratios (in this case neutron capture elements) to be then imple- 
mented in cosmological simulations (a project that we are al- 
ready pursuing). This step is crucial to avoid spending unneces- 
sary CPU time with models that do not meet the minimum re- 
quirements set by a first comparison with observations. In addi- 
tion, the comparison of our model predictions with observations 
is very useful to offer additional constraints to the theoretical nu- 
cleosynthesis predictions, which are still very uncertain for these 
elements. 

Our main results and conclusions are summarized below: 

- We first show that it is possible to reproduce the observed 
spread of [Ba/Fe] and [Sr/Fe], and the simultaneous lack of 
scatter in the a elements, in EMP stars with an inhomoge- 
neous chemical evolution model if we assume a production 
of neutron capture elements coming from: a) a relatively 
rare, but efficient, site of production (here illustrated by a 
contribution of stars in the 8-10 M© mass range), and b) a 
second site of production, which is less efficient producing 
lower amounts of neutron capture (here illustrate with a con- 
tribution from stars in the 10-40 M© mass range). More im- 
portantly, we show that despite the good agreement with the 



observed spread of [Ba/Fe] and [Sr/Fe], an extra process is 
needed to explain the observed [Sr/Ba] scatter. 

- The presence of r-process rich stars, with a common strong 
r-process signature, is an observational constraint for the first 
site of production: only if this site produces strong enhance- 
ment in relatively rare events, then the chemical evolution 
model can provide an explanation for the r-process rich stars. 
Here we show that keeping fixed this first site of production, 
but adding the contribution to the neutron capture elements 
by fast rotating massive stars (now considered as second site 
of production), then it is possible to create a scatter in the 
ratio of [Sr/Ba], as observed. 

- In particular, when considering stellar models for fast rotat- 
ing massive stars that are less conservative in their assump- 
tions (what we call here spinstar s-process models) we are 
then able, for the first time to our knowledge, to reproduce 
simultaneous the [a/Fe], [Sr/Fe], [Ba/Fe] and [Sr/Ba] trends 
and scatter observed in halo stars with an inhomogenous 
chemical evolution model. Although the proposition that the 
scatter in the observed spread between heavy and light neu- 
tron capture elements in EMP stars could be due to the con- 
tribut ion of spinstars had been made before (Chiappini et aP 
Honl), here we show for the first time quantitative estimates 
that seem to confirm this hypothesis. Notice that in this ele- 
gant solution, we are able to explain the results with just two 
different sites of production for Sr and Ba, without requiring 
more complicated scenarios. 

- The above solution is based on the existence of spinstars 
in the early Universe. Indeed, this is not the first chemical 
signature revealing the importance of fast rotating massive 
stars in the early chemical enrichment of the Galaxy. On 
the contrary, very old halo stars were previously found to 
show at least four other abundance "anomalies", which ro- 
tating models of massi ve stars can suc cessfully account for: 
rise of N/0 and C/0 (|Chiappini et al. 2006), low ^^C/^^C 
(IChiappini et al.l l2008l) and a primary-like evolution of Be 
and B (lPrantzosll2012l) . What we provide here is a 5th in- 
dependent signature, namely: an early enrichment of the 
Universe in s-process elements. 

- Our model thus predicts that stars with high [Sr/Ba] abun- 
dances should show the pollution by the s-process in spin- 
stars. Hence, these stars should present a high value of the 
even/odd isotopes of Ba (as it is typically for an s-process). 
We are aware that even with the highest resolution and sig- 
nal to noise available, the measurement of the ratio of odd to 
even isotopes of Ba could still suffer from large uncertainties 
but measurements of this ratio would be needed to confirm 
our suggestions. 

- Our results also show that to reconcile the model to the 
observations we need a r-process site of production decou- 
pled from the one of s-process. We assumed here the most 
straightforward solution, which is to consider two different 
mass ranges for the two processes, but other solutions pro- 
viding a large amount of r-process production in only a frac- 
tion (roughly 15%) of the massive stars would be equally 
valid, as in the case of ni agnetorotationally driven super- 
novae (IWinteler et al. 11201 2h . 

- At variance with suggestions so far in the literature of a 
possible "second r-process" or "weak r-process" to explain 
the observations, the production of s-process in fast rotating 
massive stars, does predict an enrichment in Ba (although not 
as strong as for Sr), and a negligible production of Eu. So, 
another natural outcome of our model is the production of a 
small spread (-0.3 dex) in the ratio of [Eu/Ba], which a sce- 
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nario involving two different r-processes cannot reproduce. 
The [Eu/Ba] spread is observed in stars with [Fe/H] 3, al- 
though more data is needed to probe this additional valuable 
observational constraint. 

We are aware that the recent simulations of star formation 
in the early univer se, find that most of the first stars are formed 
in binary systems (IStacy et al.ll2012h . However, binary stars, just 
like single stars, should also experience rotation induced mixing 
that is the key property for the production of s-process in spin- 
stars (as shown for the extreme case of long GRB progenitors by 
ICantiello et al.l (l2007l) . Indeed, the above simulations do predict 
the first stars to have large initial rotational velocities. 

As a final comment we notice that, intriguingly, the CEMP- 
no stars appear undistinguishable on the ba sis of the neutron cap- 
ture el ements from normal EMP stars. In ICescutti & Chiappinil 
our inhomogenous model with fast rotating massive stars 
was not able to account for the C/0 and N/Q ratios in the s e sam e 
CEMP-no stars. As mentioned in Cesc utti & Chiappinil (I2Q1Q). 
the explanation can reside in a pollution by a process which pro- 
duces C in a non standard chemical enrichment fashion (ie. with- 
out mixing with all the surrounding ISM). The fact that no strong 
s-process signature is observed in these stars could support the 
scenario described in IMeynet et al.l (120101) . in which a fast ro- 
tator star enrich low mass stars belonging to its forming stellar 
cluster through its stellar wind, highly enriching them in C and 
without strong dilution by the ISM. In this case, no s-process can 
pollute these low mass stars because the s-process in fast rotating 
massive stars are not removed by wind, being not present on the 
surface. However, we cannot exclude a scenario in which also 
CEMP-no are formed by an AGE companion with a dimmed 
production of s-process (although most of the CE MP-no appear 
to be single stars at the present time, see lNorris"et al 20 1 2.) . 
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Fig. 7. [Ba/Fe] vs [Fe/H] results; on the left the model for an r-process contribution only from the "standard" site, on the right the 
r-model, where both "standard" and "extended" contributions are considered. The density plot is the distribution of simulated long 
living stars for our models; the density is in a log scale, normalized to the peak of the distribution ( see bar over the Fig. [2] for the 
color scale). Superimposed, we show the abundances ratios of the halo stars (data from lFrebeill20 1 01) . The symbols are the same as 
in Fig. [3l 

We present in Fig. [71 the results of our inhomogenous model, 
with and without the "extended" site for the r-process. Thanks to 
the inhomogeneous model, it is possible to highlight that if we 
use only the "standard" r-process site (left pannel), most of the 
stars showing a low [Ba/Fe] cannot be reproduced by our model. 
This problem is not clear from the results of the homogeneous 
model because the simple line representing the homogeneous re- 
sults cannot display that a negligible amount of stars are formed 
by the model in this area. 
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